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ABSTRACT 

It has recently been argued that the observed eUipticities of galaxies may be deter- 
mined at least in part by the primordial tidal gravitational field in which the galaxy 
formed. Long-range correlations in the tidal field could thus lead to an ellipticity- 
ellipticity correlation for widely separated galaxies. We present a new model relating 
ellipticity to angular momentum, which can be calculated in linear theory. We use 
this model to calculate the angular power spectrum of intrinsic galaxy shape correla- 
tions. We show that for low redshift galaxy surveys, our model predicts that intrinsic 
correlations will dominate correlations induced by weak lensing, in good agreement 
with previous theoretical work and observations. We find that our model produces 'E- 
mode' correlations enhanced by a factor of 3.5 over i?-modes on small scales, making 
it harder to disentangle intrinsic correlations from those induced by weak gravitational 
lensing. 
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1 INTRODUCTION 



Galaxy alignments, or ellipticity correlations, have been the subject of much study over many years (Djorgovski (1987), and 
references therein), extending back well into the 1800s. Once galaxies were discovered and were found to be elliptical in shape, 
the logical next step was to see if they have a tendency to align with each other. This would provide a clue as to what they 
were and how they formed. More recently it was realised that different models of structure formation may produce different 
levels of alignment, and so searches aimed at placing constraints on these models. Searches have historically concentrated in 
nearby regions of high galaxy density (clusters and superclusters) in order to get good number statistics. The results have 
been, in gener al, inconclusive or contradic tory of previous results, although weak trends for alignments have been found for 
some samples (Cabanela & Aldering 1998). 



Recently, there has been renewed theoretical interest in long-range correlations of galaxy eUipticiti es. While initially the 
interest was in reconstructing the gravitational potential field from the orientations of galaxy spins (Lee & Pen 2000), it 
has bee n driv en lately by the increasing success of field-surveys for weak gravitational lensing by large scale structure (see 
Mellier (199!:) for a recent review). We ak lensing shear, the coherent distortion of galaxy images on the sky induced by density 



Schneider 1992 



perturbations along the line of sight (Gunn 1967; Miralda-Escude 1991; Blandford et al. 1991; Kaiser 1992; Bartclmann & 



Bartclmann fc Schneider 1999|) has now been detected by several diff e rent groups ([van Wacrbokc ct al. 200C 



Bacon, Refregier fc Ellis 2000[ [Wittman et al. 2000| ; [Kaiser, Wilson fc Luppino 2000 



Maoh et al. 2001 



Rhodes, Refregier fc 



Groth 2001). Following standard practice, all of these authors assume that all of their observed correlations in the eUipticities 
of galaxies comes from weak lensing, with an immeasurably small intrinsic signal. Intrinsic correlations, if present, may 
contaminate the weak lensing signal and so they should be considered when interpreting results from these field-lensing 
surveys. 

Theoretically, alignments are expected at some level because nearby galaxies form in similar and related gravitational 
fields; they should therefore react in similar ways to the infiuence of this field and so should have some tendency to align with 



each other. Several autho rs have revisited this question recently, with simulations (Heavens et al. 2000 ; Croft fc Metzler 200( ) 
and analytic arguments (Lee fc Pen 200C ; Catelan, Kamionkowski fc Blandford 2000; Crittenden et al. 2001). These authors 



assume that the eUipticities of the observed luminous galaxies are determined by either the shapes of their dark matter halos 
or by the halo angular momenta. In either case, the ellipticity should ultimately be determined, at least in part, by the tidal 



gravi tational field in which the galaxy forms. The existence of long-range correlations in the tidal field ( Catelan fc Porciani 
2001) will thus lead to correlations in eUipticities of widely separated sources. Observationally, Pen, Lee fc Seljak ( ^OOOD re port 



a tentative detection of galaxy alignments in the local universe using the TuUy catalog. More recently. Brown et al. (2002) find 
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a more significant correlation using data from tfie SuperCOSMOS Sky Survey^. Tfiese detections reinforce the importance of 
considering intrinsic alignments. 

An analytic theory of these correlations, even if only approximate, can provide valuable guidance on the plausible shape, 
size, redshift evolution and sample depth dependence of the effect as well as flag potential indicators of contamination. Catelan 
et al. ( 200C| ) (hereafter referred to as CKB) calculated ellipticity correlations under the assumption that a galaxy's ellipticity 
is determined by the initial halo shape. They also outlined an argument for how a galaxy's angular momentum may determine 
its ellipticity in the contex t of th e tidal torque theory using the Zel'dovich approximation. A similar model was developed more 
fully by Crittenden et al. (2001) (hereafter referred to as CNPT) who calculated correlations in ellipticity due to correlations 
in the direction of the galaxy's angular momentum vectors. Their model predicts that for low redshift galaxy surveys intrinsic 
alignments are expected to dominate over weak lensing correlations, wit h cons iderable uncertaint y in the normalisation. These 
predictions were later confirmed by the measurements of Brown et al. (2002). In a later paper (Crittenden et al. 2000) they 
also discuss ways of discriminating between intrinsic correlations and weak lensing induced correlations. 

In this paper we build upon some of the ideas presented in CKB, developing an alternate and complementary model to 
that of CNPT. We calculate in our model how correlations in ellipticity, or "intrinsic shear," can be generated by correlations 
in angular momentum, with an emphasis on seeing how large intrinsic correlations can possibly be. While some of the 
input physics of our model is the same as that of CNPT, our approach and analysis are rather different. CNPT calculate 
correlations in the direction of the angular momentum only, and for the case where the inertia tensor and tidal field tensor of 
the protogalaxy are strongly correlated. Our model assumes that thes e te nsors are uncorrelated (numerical simulations have 
recently produced evidence that this assumption is not justified, see ^2.2, although we argue that this assumption does not 
affect our results significantly), which results in a significantly different expression for the angular momentum generated. It also 
includes correlations in the magnitude of ellipticities, which may be present. The two models are normalised in different ways, 
to different physical quantities. In our calculation, we have tried to make approximations that will, if anything, overestimate 
the ellipticity correlations, in order to place upper limits on their magnitude. 

All of these models are really educated guesses as to the relationship between ellipticity and angular momentum, and are 
also chosen at least partly because of mathematical simplicity. Ideally we would like a quantitative model of galaxy formation 
which could tell us the relationship between galaxy shape/alignment and the environment in which a galaxy forms. We are 
currently far from this ideal situation, and so must make certain assumptions for our model which we believe are physically 
well motivated. It is thus both surprising and gratifying that we obtain very comparable results to CNPT. This may indicate 
that both our and their conclusions are robust to the considerable simplifications which need to be made in such calculations. 

The outline of this paper is as follows. In the next section we discuss our model for relating the angular moment um 
of a galaxy to its ellipticity in the context of tidal torque theory. In §^ we present our main calculation, explaining in §3.2 
how we normalise our model and showing in §3.3 how to relate our results to o ther measures of correlations. The reader not 



interested in the details may find the relevant equations summarized in j |3.3| . We present results for an example (ACDM) 
model in explain them in terms of the underlying physics, and consider the effects of intrinsic correlations on weak lensing 
measurements of 'cosmic shear.' Then in we discuss our results in the context of observations and previous work. We 
present our conclusions in §M. 



2 BACKGROUND 

In this section we discuss the ingredients of our model. First we introduce the concept and mathematical description of the 
ellipticity of galaxies. The key assumption we make is that the ellipticity (orientation and magnitude) is determined by the 
galaxy's angular momentum. We fix our model so that the ellipticity components transform correctly under rotations in the 
plane of the sky, and so that the magnitude of the ellipticity of galaxies with more angular momentum is larger. We implicitly 
assume that the angular momentum of the luminous galaxy is aligned with that of the dark matter halo it resides in. The 
angular momentum of a galaxy can be predicted and related to the local gravitational potential (in linear theory) using 'tidal 
torque theory'. This then allows us to calculate correlations in the ellipticity in terms of correlations in the gravitational 
potential, which is well specified in linear theory for a given cosmological model. 



2.1 Galaxy ellipticities 

It is a standard assumption in calculations of weak gravitational lensing that the source galaxies are randomly oriented with 
respect to each other. We wish to investigate to what extent this holds true. Lensing will introduce a shear in the shapes of 
the sources, so we need to describe any "intrinsic" shear that these sources have, which we do in terms of the spin-2 (complex) 
ellipticity: 

e = lele^*"* = 7i + 172 . (1) 



http:/ /www-wfau. roe. ac.uk/sss/ 
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The components 7^ are usually defined in terms of the second moments of the light distribution of a galaxy, F[x,y), as 
/ F(x,y){x^ -y^)dxdy 



71 = 



J F{x,y){x^ + y^)dxdy 



_ J F{x,y){2xy)dxdy 
J F{x,y){x^ + y^)dxdy 

It is clear that their values for a given galaxy will depend on the orientation on the sky of the coordinate axes chosen. The 
correlation functions (7^7^) are obtained by measuring 7^ at one galaxy and 7^ at another galaxy some distance away on 
the sky and then multiplying them, doing this for many galaxies separated by the same angular distance, and averaging the 
result to get the correlation function. Intrinsic shear correlations can be due to correlations in both the magnitude and the 
orientation of the ellipticity across the sky, both of which are encoded in the components (71,72)- For example, if all galaxies 
were flattened disks with axes of symmetry perfectly parallel with each other, then the correlation would be unity on all scales. 
If the disks have finite thickness, however, even if galaxies are perfectly aligned, variations in |e| with position will cause the 
correlation to vary with separation. In the opposite limit, if all galaxies are randomly oriented and have random |e|, then the 
correlations are zero and there is no alignment. 

We assume that this ellipticity is determined, at least in part, by the galaxy's angular momentum (and hence by the 
initial tidal field as we show below). For spiral galaxies this is well justified, as the disk angular momentum is perpendicular 
to the plane of the disk. There is some question as to how clo sely the angular momentum of the gas and stars in the disk 



follows that of the dark matter (Abel, Croft & Hernquist, 2001), but it seems reasonable that they should be closely related 



Elliptical galaxies tend not to have much net angular momentum, and the ellipticity seems to be determined largely by the 
velocity dispersion along the principal axes of the 3D ellipsoid. Thus our calculation may be less relevant for elliptical galaxies 
than for spirals, and the halo shapes calculation of CKB may be more applicable. 

Following the angular momentum-ellipticity discussion in CKB, using the assumption that the disk of a galaxy forms 
perpendicular to the angular momentum, the observed ellipticity of a galaxy will be 

71 = fiL,L^){Ll - lI) , 

72 = 2f{L,L.)L,Ly, (3) 

where we take the sky to be the x-y plane, and where f{L, L^) is an unknown function which determines how ellipticity scales 
with L. Since we have no firm theory to predict the form of f[L, Lz) we will take it to be a constant, C, whose value must 
be fitted empirically to the observed rms ellipticity of galaxies. As a result, in our model the ellipticity of a galaxy scales 
quadratically with its angular momentum, so galaxies with more angular momentum will be more fiattened. For the case 
where galaxies are approximated as thin disks, it can be shown that f{L, Lz) — 1/(1/^ + L1) (see CNPT), thus taking out the 
dependence on L. In this case |e| is set to a constant independent of L and only the orientation is measured by 7^. 

The price we pay for letting |e| vary with L is that galaxies in the tail of the distribution with L 3> {L) can potentially have 
an unphysical |ej > 1. We have calculated the fraction of galaxies this will apply to. In our model ellipticity is proportional to 
the square of the the gravitational potential (as we show in the next section), and so if we assume the potential fiuctuations are 
gaussian distributed, we c an c alculate the ellipticity probability distribution from this. The rms ellipticity of galaxies (which 



we take to be e = 0.4; see §3.2) gives the variance of |e|. We can thus integrate over the ellipticity probabilitiy distribution to 
find what fraction of galaxies have |e| > 1, and we find that this is ~ 4%. Larger values of |e| are exponentially rare, so they 
will not affect our calculation in a significant way. 



2.2 Tidal torque theory 



Hoyle (1949) was the first to suggest that galaxies acquire their angular momentum by tidal torques due to the surrounding 



matter distribution acting on the protogalaxy. This torque comes about from a misalignment of the p rotog alaxy's mass distri- 
bution (inertia tensor) with the local tidal field. Using the Zel'dovich approximation Doroshkevich ( 1970| ), and subsequently 
White (1984), showed that the angular momentum acquired by a protogalaxy to first order in the gravitational potential is 



iijkTjllkl, 



where Tij is the tidal field tensor at the galaxy's centre of mass and Fj is the inertia tensor of the protogalaxy. 
The inertia tensor is given by 



(4) 



(5) 



where F is the Lagrangian volume of the protogalaxy, q is the Lagrangian position vector of the constituent particles with 
centre of mass q, and po is the background density. We follow CKB in assuming each galaxy has the same eigenframe moment 
of inertia and the eigenframes are distributed isotropically. Further we suppose that two of the three principle moments 
are equal, so that the unequal moment defines a symmetry axis n = (cos 0, sin S cos 0, sin sin (^), and so the inertia tensor 
simplifies to Fj oc niUj. 

The tidal field tensor is given by 
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T,j oc aa,$(x) , (6) 

where "1> is the gravitational potential. All of the proportionality above is time-dependent, but we assume that all our galaxies 
form at a similar redshift, and so absorb all time dependent quantities into the constant C above. We normalise our theory 
in such a way that any time dependence will cancel out. 

By inspection of equation it can be seen that if I and T have the same principal axes then there is no angular 
momentum acquired to first order. There must be some misalignment between them to generate angular momentum, so the 
tensors cannot be perfectly correlated. In our calculation we make the simplest assumption, namely that they are completely 
uncorrelated. Reality lies somewhere in between the extremes of perfectly correlated and uncorrelated. Nevertheless, our 
expectation is that the large wavelength Fourier modes that are primarily responsible for the long range correlations in the 
tidal field should be statistically independent of the smaller wavelength Fourier modes that are primarily responsible for the 
inertia tensor. So the hope is that our results should not be very sensitive to this assumption. This assumption comes into play 
when we calculate the ellipticity correlation function between neighbouring galaxies. We must average over the distribution 
of inertia tensors and tidal tensors. If the two are uncorrelated we can average over inertia tensors first (as we expect the 
ellipticity correlation to be due to long range tidal field correlations), and subsequently over the tidal field. Because this 
assumption is an important part of our calculation, we discuss relevant numerical results in the next section. 



2.3 Tidal Torque Theory in Numerical Simulations 



2.3.1 Initial Angular Momentum 

Several authors have investigated the questions of a ngula r momentum generation and tidal torque theory (TTT) in numerical 
simulations. Sugerman, Summers & Kamionkowski (200C) found a relatively good correlation between the predictions of linear 
theory and the actual evolution of dark matter halos. They found that the Zel'dovich approximation we have used typically 
overpredicts the magnitude of the angular momentum generated by a factor of ~ 3, and that the spin direction is marginally 
reproduced with a larg e scatter. 

Lee & Pen (2000) differed from Sugerman et al. ( ^00C| ) in that they found quite a strong correlation between the 
predicted and simulated spin direction, again with a large scatter. They concluded that TTT is a reasonably good predictor 
of the direction of the angular momentum vector of a halo. 

Porciani, Dekel & Hoffman (2002a), in a more detailed study of TTT, used dark matter simulations to look at the spin 
(i.e. angular momentum) amplitude and direction of halos, as well as spin-spin correlations between halos. They reproduce the 
result of Sugerman et al. (2000) in finding that TTT systematically overpredicts that amplitude of the spin, but by a constant 



fraction on average, with a large scatter. The agreement between these two groups is reassuring, and gives us confidence in our 
calculation. A constant bias in the magnitude of L will be taken out by the normalisation of our model, so that the 'angular 
momentum' we calculate really does correspond (on average) to the spin amplitude of a halo. 

As regards spin direction, they find that at high redshift, TTT predictions are in good agreement with their simulations, 
although there is some scatter. They also studied spin-spin correlations of dark matter halos (i.e. correlations in the directions 
of the spin vectors among neighbouring halos) and find that at high redshift the tidal torque theory reproduces the correlations 
quite well, albeit with large scatter. These results are important for our calculation in that they show the theory on which it 
is based to be reasonably accurate, at least at high redshift. The situation is not so good at low redshift, which we return to 
below. 



2.3.2 Inertia and Tidal Tensor Correlations 



Another important effect is a correlation between the inertia and tidal field tensors. Lee & Pen (2000) found that I and 
T were quite strongly correlated in their simulations (although misaligned to a detectable degree), and this conclusion was 
supported very recently by Porciani, Dekel & Hoffman ( 2002b| ). They find that angular momentum is generated by the small 
but significant misalignment of the tensors. One way of thinking about this is to decompose the inertia tensor into a part 
perfectly correlated with the tidal tensor, and another part completely random: Itot = Icorr + Irandom (where we assume for 
now that different components of I and T are all correlated at the same level). The correlated part generates no angular 
momentum to first order so the angular momentum generated is just some fraction of what would have been generated if I 
was all random. The direction is unchanged, only the magnitude is suppressed. Such a situation does not affect our result, 
because the overall normalization of L factors out, as we show in §3.2. Catelan & Theuns ( |l99e| ) also argued analytically that 
assuming I and T are uncorrelated will r esult in an over-prediction of the magnitude of the angular momentum genera ted, 
and this has been verified in simulations (Sugerman, Summers fc Kamionkowski 2000t [Porciani, Dekel fc Hoffman 2002a| ). 

We would see an effect due to I-T correlation if the direction of L changed when these two tensors are more strongly 
correlated. One way to make this happen is if higher order terms become important when this leading order term given by 
equation (^ is suppressed. Another way is if different components of I are correlated with differing strengths to different 
components of T. Such a situation could arise if the first principal axes of I and T were strongly co rrelate d, but the second 
and third were distributed more randomly in the plane perpendicular to the first. Porciani et al. (2002b) investigated the 
correlation between the different principal axis directions of I and T. They found that each principal axis of T is strongly 
correlated with the corresponding axis of I, but all at roughly the same level (to within ~ 15%), indicating that the strength of 
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the correlation is comparable between the various components, although not identical. They also find that the initial direction 
of L (with respect to the local tidal field) for halos in their simulation is somewhat different to the direction they obtain after 
making I and T independent by randomizing the orientation of I, in that the same trends are seen but they are weaker in the 
real simulation. This could be because of the differing strengths of I-T correlations among different components, or it may be 
because higher order terms have become more important relative to the suppressed leading order term, adding larger scatter 
to the direction of L when compared to the I-T uncorrelated case. 

So where does this leave our calculation? We expect that the larger scatter and weaker correlation between L and T, 
as compared to the predictions when there is no I-T correlation, will cause galaxy alignments to be somewhat weaker than 
our prediction by some factor, making our prediction an upper limit. This is still a useful constraint, as the CNPT analysis 
explicity assumes strong I-T correlations (although see below) and so cannot be applied to the uncorrelated case to give an 
upper limit. Also this simple rescaling should be the only difference between our results due to the I-T correlation, and so 
other differences will be due to other differences in our respective models. It is important to assess which results are robust 
and which are sensitive to the simp lificati ons in both models used to calculate ellipticity. We discuss this in detail in ^ The 



differences seen by Porciani et al. ( 2002b ) are not large, however, especially when compared to the effects of late-time non- 
linear effects also seen in their simulations, so we believe that I-T correlations are not affecting the relevance our calculation 
very significantly. 

2.3.3 Late-Time Effects on Initial Conditions 

Another issue is the time dependence of a galaxy's angular momentum. The calculation presented below is a linear theory 
calculation of the angular momentum acquired in the initial stages of formation. We implicitly assume that this subsequent 
non-linear evolution of the galaxy and its interaction with its environment do not substantially alter its angular momentum. 
This is almost certainly false in dense environments such as clusters of galaxies and compact groups, where the dynamical 
time is much shorter than the Hubble time. We expect, therefore, that our analysis will not apply to them, but rather only to 
more isolated galaxies. The criterion is that they should not have exchanged significant angular momentum with neighbouring 
galax;ies/halos over a Hubble time. This is hard to quantify in practice given that galaxies form from mergers of smaller 
objects, so in a sense all galaxies have had many significant dynamical interactions in the past. Further, galaxies are not 
isolated points in the universe — they live in extended halos and it is hard to say where the halo of one galaxy stops and 
where the neighbouri ng one s start. 



Sugerman et al. (MOO) and Porciani et al. (2002a) indeed found that the spin direction of a dark matter halo changes 
with time, and that correlations between the spin direction and the initial conditions weaken significantly between the good 
agreement at z ~ 50, and z — 0. At z = there is typically a 30-40° misalignment between the initial and 2 = spin 
direction. The resulting distribution is not random by any means and clearly retains some memory of the initial conditions, 
but the spin direction does change significantly, with most of the change happening between z = 3 and z = 0. The level of 
the spin-spin correlation at z = 1 is of order 1% at a halo separation of lfe~^Mpc, which is in good agreement with work by 



Croft & Metzler (Croft & Metzler 200C) and Heavens et al. (Heavens et al. 2000[), who looked for intrinsic halo alignments in 



large N-body simulations. 

The striking feature of the Porciani et al. ( 2002a ) results is the decrease in the spin-spin correlation by a factor of 5 or 
so, as compared to the correlation at z = 50. The spin-spin correlation at 2: = 50 is in line with TTT predictions, and it is 
reall y non- linear effects such as mergers and angular momentum exchange between halos which reduce it over time. Lee & 



Pen (2000) introduced a parameter 'a' to account for the effects of I-T correlations and non- linear processes on the (linear 
theory) correlation between spin direction and the tidal field. This results in a more random spin direction, and lower spin-spin 
correlations by a factor of a^, as compared to the assumption of independent I and T and no non-linear effects. The Porciani 
et al. results indicate that it is non-linear processes acting cumulatively over time to reduce spin-spin correlations which are 
providing the dominant contribution to Lee & Pen's parameter 'a'. This strengthens our argument that I-T correlations do 
not significantly affect or invalidate our calculation, and implies that we could apply a simple correction factor to the final 
spin-spin correlation result to account for these non-linear effects. We note that it is not clear from the simulations if this 
suppression of spin-spin correlations is scale dependent or not; one might guess it is not because it seems to result from 
perturbing spin directions of galaxies in what is probably a fairly random manner. 



Also relevant is a potentially more radical idea proposed by Vitvitska et al. ( 2001 ), who investigate the possibility that 
angular momentum is generated primarily by mergers and not by tidal torques. This is perhaps not unexpected in the context 
of hierarchical models of structure formation where galaxies form by assembly of smaller parts through mergers. In this 
situation, however, the TTT can still give reasonable results even though the merger details are missed. When the tidal field 
is smoothed over a galaxy-mass region the smaller scale fluctuations that give rise to the galaxy progenitors are erased, but 
the sum of all their angular momenta (which eventually gets into the galaxy through mergers) is what is calculated by the 
theory. 

2.3.4 Summary of TTT 

There seems to be some consensus that while tidal torque theory is not perfect, it works surprisingly well and shows that 
at least the outlines of what we are observing are inherent in linear theory. That TTT works at all is somewhat surprising. 
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considering it was first thought of in the context of top-down structure formation models wherein galaxies form via the 
monolithic collapse of one large concentration of matter. For our purposes, it seems that TTT is a good starting point. To 
extend the calculation requires including non-linear effects which will add 'noise' to our predictions. These non-linear effects 
will overwrite initial correlations to some extent, and it is not yet clear whether they just add a random component to the 
angula r mo mentum or if on small scales they create their own strong alignments, or both. Indeed observations by Brown 
et al. (2002) indicate that galaxies have a fairly strong alignment in the local universe, probably significantly stronger than 
that seen between halos in simulations (at least on some scales). It is difficult to compare the strength, however, because of 
differing statistics calculated to measure correlation. 

It seems that neglecting I-T correlations affects our results somewhat, but probably not very significantly, and it is not 
clear that there is a good analytical way of taking this into account. More significant is the dilution of spin correlations over 
time which is seen in simulations. This could be dealt with by simply rescaling the analytically calculated correlations. It is 
possible (and perhaps likely), however, that since galaxies collapse to the very centres of dark matter halos, their spin direction 
is much less affected by late time effects than the outer regions of halos (where most of the angular momentum is), simply 
because there is much less of a lever arm to torque a small compact galaxy. 



3 CALCULATION 



In this section we presen t ou r main calculation. The derivation is somewhat technical and the reader interested in only the 
main results can skip to 



3.3 



Throughout we shall perform our calculation in Fourier space and express our results in terms 
of angular power spectra. This allows us to take adva ntage of much of the theo r y of spin-2 fields developed for the study o f 
polarization of e-rn radiation and gravitational waves ( Zaldarriaga fc Seljak 1997 ; Kamionkowski, Kosowsky fc Stebbins 1997 ; 



Hu fc White 1997 ). In f 3.3 we discuss how to convert our results to real space correlation functions to compare with earlier 
work. 



3.1 Power spectra 

We begin by calculating the power spectra of the ellipticities, which are each quadratic in the angular momenta (see equa- 
tion (^)). In outline, we first average our expressions for the ellipticities of a galaxy over all realisations of the inertia tensor 
to obtain an expression in terms only of derivatives of the gravitational potential. We then decompose the spin-2 ellipticity 
field into scalar (gradient or electric-type modes) and pseudo-scalar (curl or magnetic-type modes) fields in Fourier space for 
k-vectors perpendicular to the line of sight. This enables us to construct 3D power spectra for the ellipticities, which are 
convolutions over the density power spectrum. To proceed from here, we use the Limber approximation in Fourier space to 
evaluate the predicted angular power spectrum for different source distributions. It is important bear in mind that in the 
Limber approximation only modes transverse to the line of sight contribute. Evaluating correlation functions in real space is 
a relatively simple integral over the power spectra, so we calculate these for comparison with previous theoretical estimates 
and observations. 



Ignoring the normalization factors for now (§3.2), we have 



71 

72 



= 2L, 



= ninj€ijk^M ■ (J) 

As discussed above, we assume that the inertia tensor (i.e. the direction n) is independent of the tidal field. This means 
we can average over orientations of the inertia tensor separately to averaging over realisations of the tidal field. When we do 
this, we get the expectation value for a given tidal field 



(LiLj) = — {eikl^jmn^.lm^.kr. 



(8) 



This equation is the same as equation (A2) in Lee & Pen (2001), up to an overall normalisation constant. Using this, we find 



+3($^,--I'L)}, 



(9) 



71 = ~^,xx -^,yy)i^,: 

72 = {2(2$,,, 

where now 7i have been averaged over the distribution of inertia tensors. The projection of these ellipticities onto the plane 
of the sky will be a spin-2 field as in weak lensing analysis. We can construct scalar (grad or i5-mode) and pseudo-scalar (curl 
or B-mode) functions of the ellipticities by taking derivatives. We do this in Fourier space, because the analysis is somewhat 
simpler and more transparent, and in the end provides a clearer physical understanding of our results. Because the ellipticity 
is quadratic in the gravitational potential, the Fourier tran sform 7i(k) = f d'^x7i(x)e~''^'^ is a convo lution over the Fourier 



modes of the potential. The E- and _B-modes are given by ( Stebbins 199(j ; [Kamionkowski et al. 199"i ) 



6(k)r 
/3(k)fe^ 



{kx - kyh'^C^) + 2A:^fci;72(k) , 
~2kxky'yi(]s.) + (kl - fc^)7^(k) . 



(10) 
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Note that we will we working entirely in the flat sky approximation, as any correlations present will go to zero on large scales. 
So we obtain 

e(k)fc2 = ^ J ^|.(k')-l(k-k')/.(kl,kx-k'x,feL), 

= ^l>(k')4(k - k')//3(kY, kx - kl, &i) , (11) 

where k' is the convolution variable, and fe^p are functions containing all the information on the relative orientations of k 
and k' and hence the information on the derivatives of $. They are given by 

/,(a,b,c) = i(2c'-a')(6^ + (a-b)=-(axbf + 26^(a-b)) 

+^(2c^ - 6^)(a* + (a • b)^ - (ax b)^ + 2o^(a • b)) 
+3c^((a + b)^(a • b) + 2(a x b)^) , 

/^(a,b,c) = (2c'-a')(62 + (a.b))(axb) 
-(2c^ - 6^)(a^ + (a • b))(a x b) 

+3c^(a^-fe^))(axb). (12) 

Here a and b, the components of k' and (k — k') perpendicular to the line of sight, axe two dimensional vectors in the plane of 
the sky and hence their cross product is a (pseudo) scalar quantity. Because we obtained these expressions from a convolution, 
they should be symmetric under the interchange of k' and (k — k'), and this can be easily verified. It is also clear that the 
£-mode is a scalar and the B-modc is a pseudo-scalar. 

The power spectrum is defined by 

(e(k)e*(k')) = {27vf5{k - k')P«(fe), (13) 
and is given by 

(£(ki)£*(k2))A;?fci = J^(27r)^J(ki-k2)y|^P*(fc')P*(|ki-k'|) 

X /,(kl,kx-k'x,A;i) [/,(k'x,kx-kl,fc;) + /,(kx-k'x,k'x,-fci)] , (14) 

and similarly for the S-mode. Because both expressions are symmetric under the interchange of k' and (k — k') though, the 
power spectrum becomes somewhat simpler: 



(ki)£*(k2))/c?A;i = ^(277)^J(ki - k2) J ^P*(fc')P*(|ki - k'|)2/f (kl, (kx - k'x),A;i) . (15) 

We now need to evaluate the functions fcg and perform the convolution integral. Since we only consider k perpendicular 

to the line of sight, we choose our Cartesian axes such that k = fc(l,0,0), and we choose polar angles defined so that 
k' = Qfc(cos0,sin0cos</), sin^sin^) with a = k/k'. In this coordinate system, |k — k'| = + — 2ai^ (where n = cos9) 
with no dependence on 0. We find that 

[feik, a, n, (j))]^ = k^^a* [(a^ — (1 + 2a^)/Lt^ + 2an^ + (1 — + an){i — M^) cos ^^)]^ 

[fg{k,a, i^,(l>)]'^ = fc^^a'*(l — /i^)(a — /i)^(l — 2a/i)^ sin^^ . (16) 

Integrating over d(j>/{2Tr) gives functions which arc polynomials in fj, and a. Taking a factor of k^'^a!^ out of the polynomial 
and denoting the remainder by gee{cx,fi) and gi3i3{a,^i) we find 

5e.(a,M) = i((l-|-a*)(3-lV-|-19/) + 2a(3/i-2M^-17/i^) 

O 

+a{2 - 17^^ + 44/ + 19/i®) + 2a^(/u + 2/i^ - 19^®)) , 
ff/3/3(«,M) = \{a-|Jif{l-2a^lf{l-^l^). (17) 
Now the power spectrum is given by 

Pee(fe) = ^ y" ^' "2^)"'^^ -P*(Qfe)-P*(fc\/l + oP- '2aii )a'g,,{a, n) + fc,, -terms. (18) 

where A:||-terms indicates terms with line-of-sight component 7^ which will not contribute to our final result. The cross 
power spectrum Pepik) is identically zero as it is parity violating. We use Poisson's equation, and express the ellipticity power 
spectrum in dimensionless form as 
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and similarly for the B-mode power spectrum. Here is the contribution to the mass variance per log k 
da^ _ k^P(k) 



dlog k 



(19) 



(20) 



To proceed to the angular power spectrum we use the Limber approximation. In Fourier space the Limber approximation 
states that the angular power spectrum of a projection, J dx w(x)S, of a scalar field ^(x), is (Kaiser 1992; White & Hu 200C) 



e{2e + i)Ce _ TT 

4^ ^ I 



Xdx w (x)Agg(fcx = ^,a) : 



(21) 



where w(x) is the weight function at radial distance x- Before we compute our final result we need to return to the issue of 
the normalization constant C. 



3.2 Normalisation 

We determine C empirically by fitting the model prediction for the mean-square source ellipticity to that which is observed, 
= (7? + 72) 

= C^iiLl + Llf). (22) 

We are assuming here that all of a galaxy's ellipticity is due to its angular momentum acquired via tidal torques. It is likely 
that some fraction of the ellipticity is due to other effects (e.g. halo shapes as in CKB) in which case our normalisation 
provides an upper limit to the angular momentum correlations. 
Averaging over all orientations of L we find 

e- = ^(L^> = ^(Ly. (23) 



The expectation value of L is calculated (Catelan & Theuns 1996) to obtain 



where all the time dependent terms and the magnitude of the inertia tensor have been absorbed into the constant C. Note 
C cx 1/(L^) so any time dependence in L cancels out in the expressions for 7^. We can therefore calculate everything according 
to linear theory extrapolated to the present {a — 1). Also, the magnitude of L is irrelevant in our calculation; what is important 
is the ratio of to the mean squared value {L^)- The quantity, 

/f A?,,.,.. 1)(JMM)'. (2.) 

is the variance of the density field smoothed on a scale R, associated with the size of a galaxy mass regi on in the initial density 



field. We will take R = Ih ^Mpc , but calculate the dependence of our results on this choice in §4.1. Finally we obtain the 
normalised ellipticities, 



71 



1125 ijLl-Ll) 
32 {In^HifaHR) 



,1125 e2LxLy , , 

We take e — 0.4 to present our results, as it is a good fit to the observed rms galaxy ellipticity (see CNPT). 



3.3 Summary 



We are now in a position to derive our central result. Including the normalisation factors from §3.2, we obtain the angular 
power spectrum 



^(2^ + l)C7r nf f da , f\ 

TT / Xdx-U} (X) / — A„( — ) / dfi — g,,{a,ii), (27) 



47r 32a^{R) £ J ^ 'J a x 7_i (1 + - 2a^i)7/ 

and similarly for the power spectrum. The cross-spectrum vanishes due to parity. 

In integrating equation ( p^ ) the factor (1 + — 2afi)~'^^^ must be treated carefully near 2> 1 and k ~ k' . We make an 
approximation in this region which is good to < 7 per cent in the worst case (a low redshift galaxy distribution and £ > 10^). 
The weight function depends on the number of sources per unit redshift. We take 



© 0000 RAS, MNRAS 000, 000-000 



Theoretical Estimates of Intrinsic Galaxy Alignment 9 
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10 100 1000 10* 10 100 1000 10* 10 100 1000 10* 
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Figure 1. The various angular power spectra discussed in the text for our fiat ACDM model. The solid line is the predicted weak lensing 
signal , the dashed line the "intrinsic" ellipticity EE power spectrum and the dot-dashed line the intrinsic BB power spectrum. The 
first panel is for a low redshift source galaxy distribution of mean redshift {z) = 0.1 using the redshift distribution described in the text. 
The second and third panels are for (z) = 0.3 and 1.0 respectively. 



dn 



oc X exp 



(28) 



with J dn — 1. The case a — 1, (3 = 4 approximates the source distribution of a flux limited survey, and we adjust x* to 
obtain a desired mean source redshift (see below). 

It is straightforward to go from the angular power spectrum to the correlation function. The two are related by a slight 
generalization of the Hankel transform as 



<Ti7i) = / -^(Cr [Jom + M£9) COS' 



+ Cf [Jo{£e) - Mie) COS 40]) , 



(7272) 



ld£ 
4% 



(Cf [Jo{£e) ~ J4{£e) cos4<^] + Cf*^ [Jo{£0) + JaW COS 40]) , 



(29) 



where the separation between the two points has polar coordinates {9, 0). For comparison with observations we also calculate 
the "ellipticity variance" smoothed on some angular scale 9. This is obtained by convolving the ellipticities 7i with a filter 
function and taking the autocorrelation of this at zero lag. For a top-hat real space circular filter of radius 9, the variance is 
given by 



a\9)^ / ^-§{CT+C, 



00\ 



2Ji( 



£9 



(30) 



4 RESULTS 

To demonstrate our results we use a flat ACDM cosmological model with para mete rs ($7™. flb/i^, f^A, ft, n, erg) = 
(0.3,0.018,0.7,0.67, 1.0,0.9), with the fit to the transfer function of Eisenstein & Hu ( [l997| ). Fig. |l| shows the ellipticity- 
ellipticity power spectra, along with the power spectrum predicted by weak lensing for three source distributions (with mean 
source redshifts (^arc) = 0.1, 0.3, and 1.0). 



4.1 Amplitude of the angular power spectra 

We first discuss the overall amplitude of the intrinsic and lensing power spectra. The scales of most relevance for recent 



measurements of intrinsic correlations (Brown et al. 2002) and for current weak lensing surveys are from £ ~ 100 (S ~ 2°) to 
£ ~ 1000 (9 ~ 10'). For our low redshift galaxy sample with (^src) = 0.1 the intrinsic signal dominates the weak lensing signal, 
while the situation is reversed for the higher redshift samples. This is one of our main results, and agrees with calculations 
and measurements by previous authors (CKB, CNPT, Heavens et al. 2000, Croft & Metzler 2000, Brown et al. 2000). Our 
result adds further support that this is the order of magnitude to expect intrinsic correlations due to correlations in the tidal 
field. 



The result can be easily understood as the combination of two effects (these were also alluded to by Miralda-Escude (1991) 



as an argument for ignoring intrinsic alignments in deep surveys). For lensing the power spectrum is proportional to the square 
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of the projected mass density, which increases with increasing survey depth. Intrinsic correlations, however, are located in the 
source plane due to galaxies physically close to each other. For deeper surveys, the probability of two galaxies close to each 
other on the sky being physically close to each other is much smaller than for shallow surveys. So the intrinsic correlations 
get 'washed out' by this projection effect. 

This result shows that for deep weak lensing surveys with (zsrc) ~ 1.0, intrinsic correlations are unlikely to contribute to 
the signal at a significant level, as long as the redshift distribution of sources is broad. For shallower surveys, however, they are 
expected to be the dominant signal. This has been convincingly demonstrated by Brown et al. (2002) using the SuperCOSMOS 
Sky Survey. They detect a significant intrinsic correlation of the order of magnitude of our results. We compare our results 
with these observations in more detail in the next section. 

There are two main uncertainties in the amplitude of our prediction. The first is that for very large £, we expect non-linear 
clustering of galaxies to erase the initial correlations due to interactions with neighbouring galaxies. This may randomise 
ellipticities altogether, or it may create new correlations with a different amplitude and scale dependence. Also, based on 
numerical results, the overall amplitude of the signal may be scaled down by some factor, due to the weakening of spin-spin 
correlations over time. Secondly, the smoothing scale in our normalisations affects our results. We use R — l^~^Mpc as the 
comoving scale containing on average a galaxy mass (W^^Mq) at the background density. Our results scale as a~'^{R), so 
changing the smoothing scale alters things significantly. For our cosmological model, with as = 0.90, we find that (T(i?) — 3.6, 
2.8 and 2.0 for R — 0.5, 1.0 and 2.0/i^^Mpc respectively. With these numbers, increasing (decreasing) our smoothing scale by 
a factor of two increases (decreases) our results by a factor of 3.5 (2.9). 

The rms ellipticity of galaxies empirically fixes the normalisation. We choose e = 0.4 as it is unlikely to be larger than 
this, but it may be somewhat smaller. As stated earlier, it is also possible that only a fraction, / < 1, of the rms ellipticity is 
determined by the angular momentum, with the rest contributed by other processes (e.g. halo shapes). In this case our results 
should be multiplied by this factor of /, thus lowering the predicted power spectrum. We show all our results assuming / = 1 
to emphasize how large intrinsic correlations could be. 



4.2 Power Spectrum Shape 

The shape of the intrinsic power spectra is similar for all three redshift distributions, and is refiective of the underlying 
physics. On large scales (small £) the (log) slope is 2 for both the EE and BB power spectra, consistent with shot noise 
(C" — Cf^ ^constant). This is expected as there should be no intrinsic correlation between sources that are widely separated 
on the sky. The lensing signal falls off more slowly at large angles because it is an integral over the power spectrum whereas the 
intrinsic signal is an integral over the square of the power spectrum. Also lensing occurs between us and the source galaxies, 
whereas intrinsic correlations are at the source. 

The power begins to fall significantly below the shot noise power on an angular scale corresponding roughly to the turnover 
in the mass power spectrum at the mean source redshift. The linear theory power spectrum changes slope at ~ 0.02/i Mpc~^, 
and fcx ~ ^ at the mean redshift, so we can work out what value of I we expect this feature in A^(fc) to show up. Mean 
redshifts oi z — 0.1, 0.3 and 1.0 correspond to distances of x = 280, 730 and 1650/i~^Mpc in our cosmological model. Thus 
we expect the roll over to occur at roughly ^ ~ 6, 15 and 35 for the three source distributions in Fig. |l| This is indeed what 
we see in that the E- and B-modes start to differ from each other and to change slope for larger values of I. 

The large I slope flattens off (and may even start to decrease for the low redshift population) due to the fact that the 
3D matter power spectrum becomes flat at large k. We emphasize again, however, that on very small angular scales our 
predictions may be completely overwritten by non-linear clustering effects. 



4.3 Ratio of E- to B-mode power 

The i5-mode intrinsic signal is enhanced over the _B-mode by a factor of C"/Cf ~ 3.5 on small scales (^ ^ 1). This goes 
to a constant ratio for all three galaxy redshift distributions, just shifted to smaller angular scales for deeper surveys. On 
very small scales the power spectrum becomes roughly constant (increases only as logfc), so this ratio must be determined 
by the geometry of the problem and the tensor nature of the source. A similar result was obtained by (Hu & White 1997) 
in studying the polarisation induced in the CMB by tensor perturbations (gravitational waves), where they found that the 
ratio tended to C"/C^'' = 13/8 in the small angle limit. Given the mathematical similarity between polarization and galaxy 
ellipticities, it is likely that similar geometrical considerations are producing the B-mode suppression in our model, although 
the suppression is stronger in our case. 

An intuitive way of thinking about this B-mode suppression is that an isolated point mass can generate only B-modes. 
We believe that this is the reason our B-modes are suppressed on smaller scales where the density held can be increasingly 
described in terms of "objects". As shown by Lee & Pen (2001), if the tidal field has ordered eigenvalues Ai > A2 > A3 along 
its principal axes then (in this principal axes frame) Li oc (A2 — A3)/23 and similarly for the other components. If I and T 
are uncorrelated we expect L2 to by largest from the ratio of the sizes of the eigenvalues. Porciani, Dekel & Hoffman (2002h) 
showed that in their simulations, in general L\ is smaller than the other components, so the angular momentum is largely in 
the plane perpendicular to the first principal axis of T. This can be understood by considering a proto-galaxy forming in the 
vicinity of a large point mass (or equivalently a large spherically symmetric mass distribution). In the limit that the tidal field 
is dominated by this mass, the flrst principal axis of T will lie in the radial direction (i.e. along the separation vector between 
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Figure 2. Correlation functions and comparison with other work, all for the low redshift (z) = 0.1 source distribution. Note that 
different quantities are plotted on each panel, with a different y-axis scale on the first one. Left panel: Comparison of the correlation 
functions for intrinsic ellipticity correlations (in our model) and lensing induced correlations, as calculated from the power spectra in 
Fig. 1. Middle panel: Comparison with the theoretical results of CNPT. Plotted is the sum of the two correlation functions in each case. 
Curve is from Fig. 5 of CNPT for a = 0.24. Note the similarity of our results, although the amplitude of our result is slightly larger. 
Right panel: Comparison with the observations of Brown et al. (2000). Plotted is the ellipticity variance as described in the text. The 
vertical line at 10' corresponds to a linear scale of roughly l/i~^Mpc at the mean source redshift, below which we do not expect our 
calculation to be relevant. 



the proto-galaxy and the point mass), and the second and third will be in the tangential directions. The second and third 
eigenvalues are identical, because there is no preferred tangential direction, making Li = 0. Thus the angular momentum 
acquired due to tidal torques will be perpendicular to its separation vector from the point mass, and hence the galaxy will be 
oriented radially. This is a pure i?-mode pattern. To see this in the context of our calculation, notice that only the radial first 
derivative of $ is non-zero, so all cross derivatives (with i 7^ j) vanish. These cross derivatives become cross products of 
fc-modes in Fourier space, so inspection of equation ( |l2| ) shows that the B-mode must vanish. So if the shear field is dominated 
by isolated point masses, then the correlations will be mostly i5-modes. 

That our model produces B-modes at all is due to the ellipticity components 7^ being quadratic in L and hence in $(x). 
This can be seen in our expression for /3(k) in equations ( |ll| ) and ( [I4 ) as a convolution over the Fourier modes of the potential. 
Examination of equations ( pT|Jl^ ) shows that most of the contribution comes when k and k' are quite different in magnitude. 
We interpret this as indicating that large-scale perturbations on the small-scale potential field generate the B-modes. 

One of the main ways of discriminating between weak lensing and systematic effects is that weak lensing produces only 
iJ-mode power. Shot noise has C" = = constant, so looking at the B-modes gives an estimate of the level of noise in 
the iJ-mode lensing signal. Our model indicates that there can be significantly larger intrinsic i5-mode power than would be 
naively be expected from looking at the level of B-mode power. Thus the suppression of the B-modes in our model means 
that intrinsic correlations may be hidden in a low signal-to-noise detection of £-mode power. On the other hand, this factor 
of 3.5 suppression should be an upper limit to a real measurement. Adding shot noise, with C" = Cf i will only reduce the 
E- to B- ratio. Non-linear effects adding small-scale correlations may also alter the ratio. 



5 DISCUSSION AND COMPARISON WITH OTHER WORK 

5.1 Correlation Functions and Comparison with CNPT and Simulations 

We calculated correlation functions from the power spectra in Fig. |^ using equation (|2^) for the low redshift source distribution. 
These are plotted in the left panel of Fig. ^. In evaluating equation ( p9| ) we use = 0, corresponding to the separation vector 
of the two galaxies being along the a::-axis. Note that the {7272) function goes anticorrelated for large separations, and so we 
plot its absolute value. The physics producing the curves is less clear than in the power spectrum analysis, but they do show 
the same basic features — intrinsic correlations should dominate weak lensing signal at low redshift. 

Computing correlation functions allows us to compare our results with the work of CNPT who pursued a real space 
correlation function analysis of correlations in the orientations of galaxies. This comparison is shown in the middle panel of 
Fig. where we are comparing with the a = 0.24 curve in Fig. 5 of CNPT. Q The first thing to note is that our results 
are broadly similar to theirs. The amplitude of our signal is slightly higher, although increasing the value of the correlation 



Note that CNPT seem to define 'a' differently to Lee & Pen (2000), so that the a = 0.24 they plot (and which we compare to) may 



t 

not correspond exactly to the a = 0.24 found by Lee & Pen in their simulations. 
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parameter 'a' in their model can reverse this. This shows that extrapolating their model to a = 1 (for the case of uncorrelated 
I and T tensors) overestimates the correlations. The reason for this is an approximation made in their calculation which 
assumes a ^ 1, so it was not designed to work when a — 1. It is interesting that our results only differ by ~ 50 per cent 
at ^ ~ 1'. We would expect, because they normalise their model to take account of the suppression of spin-spin correlations 
seen between dark halos in numerical simulations, that their result should be significantly smaller than ours. It may be that 
differences in the normalisation of the underlying dark matter power spectrum, or other model differences, are decreasing the 
difference between our results. Overall, given the major simplifications and assumptions that go into each model, we feel it is 
encouraging that both models obtain similar results. 

The slopes of our correlation functions are also somewhat different, with our correlation function falling off more slowly at 
large angular separations. We suspect that this is probably reflecting the different matter power spectra used in the calculations 
(CNPT used a power law whereas we use the full linear-theory power spectrum appropriate for our cosmological model), but 
it may also be due to differences in our models. The fact that our results are so similar adds support to the calculation of 
CNPT, and suggests that both their results and ours are reflecting the physics of tidal torques and not the details of our 
respective models. 

This conclusion is further supported by recent numerical work. Heavens et al. ( ^000| ) assumed a galaxy disk will form 
perpendicular to the angular momentum vector of the halo it lives in and calculated spin-spin correlations from large numerical 
simulations. They found significa nt sp in correlations for halos at z = 1, and less significant correlations at low redshift. At 
the same time. Croft & Metzler (2000) measured halo shape correlations in similar large sir nulati ons and found correlations 
of similar magnitude to (but slightly stronger than) the spin correlations of Heavens et al. (2000). Our results are in rough 
agreement with the amplitude of the signal found in these numerical simulations. 



5.2 E- and B-Modes 



In a second paper, Crittenden et al. (2000) calculate E- and _B-mode correlation functions by constructing local estimators 



of the E- and _B-mode contributions to the ellipticities in real space. While their method of analysis is very different to the 



power spectrum approach we have presented here (see Hu & White (2000) for an application of this method to simulated 
weak lensing data), the calculations are mathematically equivalent. We find that our results are qualitatively different from 
each other. For two different matter correlation functions (power laws with slopes —1 and —3/2) they find that the E- and 
-B-mode correlation functions are the same for small separations in both cases. This is to be contrasted with the power spectra 
we obtain in Fig. |l|, where the two are the same on large scales and differ on small scales. They also find that on large scales, 
for ^(r) oc r~^, the i?-modes are enhanced somewhat over B-modes, while they are identical for ^(r) oc r"''^^. The fact that 
the relative contributions change depending on the power law slope of the correlation function indicates that this is the source 
of the difference between their results and ours. We use a power spectrum which has different large and small scale slopes to 
the corresponding slopes in their correlation functions, so our resulting E- and B-modes look different in the large and small 
angle limits. We have investigated the power spectrum dependence of the E/B ratio by considering 'tilted' power spectra with 
scalar perturbation spectral indices of n = (0.75, 1.25) as well as the canonical n = 1 which we used to present our results. We 
find that indeed the ratio (at £ — 5000 for (2) — 0.1) decreases with increasing spectral index. A 25% change in n produces 



roughly a 12% change in the E/B ratio. The reason for this is alluded to in §4.5 , where we showed that for a given value of £, 
the contributions to the E- and B-modes are from a range of different scales, and the ranges are different for the two modes. 
Tilting the power spectrum changes the ratios of power on different scales, and so will change the relative contributions to 
the E- and _B-modes. 



5.3 Comparison with Observations 

In the right panel of Fig. H we p lot the intrinsic ellipticity variance as calculated from equation (^) together with the 
observations of Brown et al. ( [200!4 )- Observations by Pen et al. ( ^000| ) are at a similar level, but at lower statistical significance 
and over a smaller range of scales, so we only compare to Brown et al. here. At first sight it appears that our calculation does 
not fit the data very well. It is encouraging, however, that we have the right order of magnitude, and there are a number of 
other factors to take into account. At redshift z = 0.1, 10 arcminutes (indicated by the vertical line in the figure) projects 
onto a linear separation of ~ 1/i^^Mpc. We expect that non-linear clustering and dynamical interactions will have erased 
initial conditions for galaxy separations much smaller than this. Non-linear effects may produce their own alignments, and we 
suspect this is what is giving the large observed signal at small separations. 

On larger separation scales, we have argued that our prediction should be an upper limit to correlations due to tidal 
torques. So the fact that we have obtained the right order of magnitude and are on the high side of the data is a good sign. 
The agreement of our results and the data on angular scales from 10' to 100' is really quite good, taking this into account. 
In order to conclusively test whether this observed correlation is due to correlations in the tidal field we need data that goes 
out to larger scales. It is important to see whether the variance stays constant as the data hint at, or whether it falls off on 
larger scales as the theory predicts. 
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6 CONCLUSIONS 

In conclusion, we have presented a new analytic calculation of ellipticity correlations due to correlations in the tidal field 



galaxies form in. This correlation is probably giving the observed signal seen by ( prown et al. 2002[ ) on scales of 10' to lOO'. 
Our calculation is complementary to previous work of CNPT, who performed a similar calculation with similar input physics 
but difi'erent assumptions relating the physics to the observed ellipticity. The similarity of our results adds support to their 
validity, and indicates that they may be reflecting the physics of tidal torques rather than the detailed model assumptions. 
We extend previous work in applying the Fourier space E- and B-mode power spectrum formalism to the problem (but 



see Crittenden et al. (2000) for an application of this decomposition in real space). We find that i?-modes are enhanced by a 
factor of ~ 3.5 over B-modes on small scales in our model, which can be understood intuitively by noting that isolated point 
masses can generate only _E-modes. This means intrinsic i?-mode contamination of weak lensing signal can be c onsid erably 



larger than that implied by looking at the B-modes. This result differs from the findings of Crittenden et al. (2000), who 
found that E- and B-modes have the same amplitude on small scales. This is probably mostly refiecting the different matter 
power spectra used in the calculations, but may also be at least partly due to differing model assumptions. 

We confirm the results of previous authors in finding that intrinsic correlations are a small but possibly significant 
contaminant for weak lensing surveys with mean redshifts of order unity. We have shown that the contamination can be as 
large as 30% (in rms correlation) on degree to 10-arcmin angular scales, but this should be an upper limit according to our 
calculation. Thus, while at present uncertainties due to intrinsic shape correlations are at most comparable to (and probably 
smaller than) other observational uncertainties for deep surveys, they will become important for future precsion surveys 
designed to put strong constraints on cosmological parameters. They will also be important for s urveys in which the galaxy 



sample is distributed into narrow redshift bins by photometric techniques. Croft & Metzler (200C) have shown that reducing 
the width of the redshift distribution increases the intrinsic signal while leaving the lensing signal essentially unchanged. 

Shallow surveys, on the other hand, are a very good probe of intrinsic correlations. These correlations should be investi- 
gated in more detail because they are interesting in their own right. Their existence and strength may be able put constraints 
on models of galaxy formation, and on the relationship between baryonic matter and dark matter halos by further comparison 
with numerical simulations and analytic theory. 
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